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ELECTRON DYNAMICS IN A PLASMA FOCUS 


Frank Hohl , S. Peter Gary,* and 
Patricia A- Winters 
Langley Research Center 


SUMMARY 

Results are presented of a numerical integration of the three-dimensional 
relativistic equations of motion of electrons subject to given electric and mag- 
netic fields deduced from experiments. Fields due to two different models are 
investigated. For the first model, the fields are those due to a circular dis- 
tribution of axial current filaments. As the current filaments coliapse toward 
the axis, large azimuthal magnetic and axial electric fields are induced. These 
fields effectively heat the electrons to a temperature of approximately 8 keV 
and accelerate electrons within the radius of the filaments to high axial veloc- 
ities. Similar results are obtained for the current-reduction phase of focus 
formation. For the second model, the fields are those due to a uniform current 
distribution. Both the current-reduction and the compression phases were stud- 
ied . There is little heating or acceleration of electrons during the compres- 
sion phase because the electrons are tied to the magnetic field. However, dur- 
ing the current-reduction phase, electrons near the axis are accelerated toward 
the center electrode and reach energies of 100 keV. The computations are in gen- 
eral agreement with experimental results. A criterion is obtained which limits 
the runaway electron current to about 400 A. The criterion depends only on the 
electron temperature and thus is subject to experimental observation. 


INTRODUCTION 

Plasma-focus devices have been extensively investigated in a number of lab- 
oratories. (See refs. 1 to 8.) Figure 1 illustrates the coaxial plasma-focus 
apparatus. The region of interest in the present investigation is the small 
cylindrical shaded region directly over the center electrode. This is the 
region where the dense, high- temperature plasma appears in the final stages of 
focus formation. Operation of the plasma focus with deuterium as the filling 
gas yields intense bursts of neutrons and X-rays. A two-dimensional fluid model 
of the plasma focus has been developed by Potter (ref. 9), and the model 
described well the initial phases of focus formation. However, such fluid mod- 
els do not include the detailed processes of ion and electron heating occurring 
during the final phase of focus formation. To investigate the process of ion 
heating and acceleration in the plasma focus, a number of models have been used 
to determine ion trajectories in given fields deduced from experimental data. 
(See refs. 10 to 12.) Some of the more recent work on ion dynamics based on an 
improved model (refs. 13 and 14) is in good agreement with the experimental 
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data. The purpose of the present paper is to present results on the electron 
acceleration in the plasma focus and to discuss the associated X-ray production. 

X-ray emission from the plasma focus has been investigated by many differ- 
ent experimental groups. (See refs. 3, 4, and 15 to 19.) The experimental 
results show that the plasma focus emits a substantial flux of X-rays with ener- 
gies greater than 100 keV. Also, the hard X-ray intensity is greatest in the 
direction pointing from the plasma to the anode, indicating that the X-rays are 
produced by high-energy electrons striking the anode (ref. 19). 

Newman and Petrosian (ref. 20) have proposed a model of the plasma focus 
where electrons are accelerated during the current-sheet collapse (compression) 
phase, which would accelerate electrons away from the anode to produce X-rays in 
the plasma. However, this conclusion is not supported by experimental evidence. 
For example, figure 2 shows a picture from an X-ray pinhole camera of the plasma- 
focus formation process from the 50 kV 50 kJ plasma focus at the Langley Research 
Center. The usual pinhole camera picture shows only the X-ray emission from the 
electrode and from the very dense cone-shaped plasma directly above the center 
of the electrode. This result is caused by electrons accelerated downward and 
striking the electrode or the dense plasma. To determine the presence of any 
upward accelerated electrons, an aluminum plate was placed above the center elec- 
trode as shown in figure 2. No X-ray image of the aluminum target was observed, 
thus indicating that electrons are accelerated primarily downward toward the cen- 
ter electrode and not upward . 

The uniform-current model used in the present paper shows that electrons 
are accelerated toward the center electrode in agreement with experimental 
results. The plasma-focus parameters used in the present paper are the same as 
those used in an earlier paper on ion heating in a plasma focus (ref. 14). That 
is, the peak current in the focus is taken to be 1 MA. During the collapse 
phase the collapse velocity is 180 km/s and the current-reduction phase results 
in a 30-percent decrease of the current in a time of 50 ns. 

Two models of the current distribution are studied in the present paper. 

The first distribution is that due to a circular distribution of current fila- 
ments as observed by Bostick et al . (ref. 21). Evidence of such filaments has 
also been observed by radial tracks on the copper electrode from the plasma 
focus at the Langley Research Center. An example of such tracks is shown in fig- 
ure 3- The central hole in the electrode is the result of energetic particle 
interaction resulting from focus formation. The second model investigated is 
a uniform current distribution used previously for the ion dynamics studies 
(refs. 13 and 14). For both models the collapse and the current-reduction 
phases are studied . 


SYMBOLS 

A magnetic vector potential, webers/meter 

B magnetic field, teslas 

c speed of light, meters/second 
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r, e,z 


electric field, volts/meter 

magnitude of electronic charge, coulombs 

current, amperes 

runaway current, amperes 

current density, amperes/meter^ 

mass, kilograms 

electron mass, kilograms 

ion mass, kilograms 

number of current filaments 

electron-number density, meters"^ 

radius of current filaments or radius of current-carrying region, 
meters 

initial radius of plasma column, meters 
cylindrical coordinates 
position vector 

radius of outer electrode, meters 
runaway radius, meters 

temperature derived from -1 m^(V - <V>) 2 ), electronvolts 


effective energy of electron beam derived from 
electronvolts 



time , seconds 
velocity 

velocity vector, meters/second 
average drift velocity, meters/second 
collapse velocity, 180 kilometers/second 
thermal velocity, meters/second 
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x,y,z rectangular coordinates 

3 - y Q en e V av /2 

\x Q permeability of free space 

x time constant per second 

cyclotron frequency per second 
Subscripts : 
e electron 

i ion 

max maximum 

r,0,z r- , e-, and z-components 

x,y,z x-, y-, and z-components 

Notation : 

< > root-mean-square value 

^ unit vector 


THE MODEL 


Equations of Motion 


The 

magnetic 


relativistic equations of motion for electrons subject to electric and 
fields are 


and 



(D 


^ = V(t) (2) 

dt 

By neglecting the displacement current, the relevant Maxwell equations are 

V x B = Mo J (3) 
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and 


V x I = - -4 (4) 

9t 

In cylindrical coordinates with the current in the z-direction, the fields 


B = 0B 0 + rB r 


E = zE z (6) 

Equations (1) and (2) are integrated numerically by using the time-.centered leap- 
frog method . 

Current Filament 

The vector potential for a ring of current filaments in the plasma focus is 
approximated by the vector potential for an isolated ring of n current fila- 
ments equally spaced around a circle of radius R. That is, 


u I / 
A-(r,e) a - -2- in { 


R 2n _ 2R n r n cos n9 + r 2n ) + C 


where I is the current^in each^filament and C is a constant. The magnetic 
field is obtained from B = V x A, or 


B e = 

0 3r 


(S) °°° - - (S)‘ 


2 irr /v >\ n / w, v ^ 

1 - 2(E) cos ne * (E) 


B_ . 1 . 


The magnitude of B is 


y Q nI y r ) sln n9 

2ltr - . ,»>.K /r\ 2n 

1 - 2(E) oos ne * (E) 


B = Wb 2 ,B p 2.V£ 
V e r 2irr 


1 - 2(±) cos 


- 9 * (fy 


For the parameters n = 10, R = 10“3 m, and nl = 10^ A, the constant B 
contours are shown in figure 4. There are large gradients of B near the 
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filaments and B is essentially zero in the region interior to the ring of 
filaments. This behavior is again displayed in figure 5 which shows B as a 
function of r. The electric field induced by the changing current and resul- 
tant magnetic field is obtained from 



During the compression phase the radius R of the current filaments is taken to 
collapse with the experimentally determined velocity V c = -dR/dt of 180 km/s. 
(See refs. 4 and 14.) Thus, the induced electric field during the compression 
phase is 


E z = 


, Vi 

C 2irR 


1 


/ r \ 

1 - (-J cos ne 

2 (f)” 008 88 ♦ (sf 


( 12 ) 


Note that for r = r Q (where r Q = 0.05 m, the radius of the outer electrode), 
E z is essentially zero and no term need be added to equation (12) to satisfy 
the boundary condition E (r ) = 0. Figure 6 shows the radial variation of the 
induced field at various angles between two filaments. Note that the angle 
between filaments is 36° and the angles 0° and 36° pass through the center of 
filaments . 


During the current-reduction phase, R remains constant and the electric 
current I varies as 

I h I Q ( 1 - xt) ( 13 ) 

where Iq = I (t = 0) and where t is determined in agreement with experimen- 
tal results to obtain a 30-percent current reduction in 50 ns. Thus, the axial 
electric field during the current-reduction phase is 

E z = - - — [in (R 2n - 2R n r n cos n6 + r 2n ) - In (R 2n - 2R n r Q cos n0 + r Q 2n )j 

( 14 ) 

The second term in the equation is a constant such that the induced elec- 
tric field becomes zero at radius r Q which is the radius of the outer elec- 
trode. Figure 7 shows the radial variation of the induced axial field during 
the current-reduction phase as given by equation (14). Note that the 50-ns 
current-reduction phase is probably too long. Thus, the induced field in the 
central region during the current-reduction phase is expected to be about the 
same magnitude as the peak field during the compression phase as shown in 
figure 6. 


Uniform Current 


The magnetic and 
in the present report 


electric fields for the uniform-current density model used 
are the same as those used previously in reference 14, and 
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their derivation will not be repeated here. During the compression phase the 
magnetic field is 


B (r,t) = 


_ 


2irR‘ 


(r < R(t) ) 


B 0 (r,t) = 


2irr 


(r > R(t) ) 


where R(t) is now the radius of the current-carrying region and I is the 
total current. The corresponding induced electric field is 


E (r,t) = -^^(r 2 - R 2 ( t ) ) 


2ttR' 


E (r,t) = 0 


(r ^ R(t) ) 
(r > R(t) ) 


For the current-reduction phase the magnetic field is 


B 0 ( n ,t )' 


t J o I 0 

2irR 2 


r( 1 - xt) 


(r ^ R ( t ) ) 


B ( r ,t ) = ^_0(1 - xt) 
0 2irr 


(R(t) < r < r 0 ) 


resulting in an induced electric field given by 


E ( r) = - ^^fr 2 - 2R 2 


4ttR c 


In 


(r < R) 


E z (r) 


-^(l + 


2 In r - 2 In - Q 


(R < r < r Q ) 


RESULTS AND ANALYSIS 
Current-Filament Model 

Filamentary current structure has been observed during the collapse or com- 
pression phase of focus formation. It, therefore, seems appropriate to investi- 
gate the effect of the associated fields on the plasma electrons. Details of 
the fields are shown in figures 4 to 7 . Figure 8 shows the evolution of the 
electron distribution represented by 3000 electrons during the compression 
phase . The initial conditions of the electrons are generated by a pseudo random 
number generator. The initial velocities were given a Gaussian distribution, 
whereas the initial positions were given a uniform distribution. The magnetic 
and induced electric fields acting on the electrons are given by equations (8), 
(9), and (10). The variable R defining the radial position of the current fil 
aments is given by 
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R < t ) n R 0 - V c t 

Experimental results (ref. 4) indicate that the aforementioned expression is 
valid for 0 ^ t i 50 ns. However, because of the rapid electron response to 
the strong induced electric fields, the calculations are performed for only a 
fraction of a ns. Note that electrons are expelled from the regions of high 
magnetic field near the filaments. Also, electrons concentrate in magnetic 
wells between filaments and in the central region. The evolution of the total 
electron velocity as a function of radius is given in figure 9. As can be seen 
from figures 5 and 6, the magnetic field is essentially zero in the region inte- 
rior to the ring of filaments, whereas the axial electric field is large in this 
region. Thus, the electrons in the central region quickly reach large axial 
velocities and essentially free stream under the influence of the induced axial 
electric field. As the electrons reach an axial distance of 1.5 cm from their 
initial z = 0 position, they are assumed to have escaped the focus region and 
are taken out of the calculations. This explains the absence of accelerated 
electrons in the central regions after t = 0.15 ns in figure 9- Figure 10 dis- 
plays the energy of the electrons during the compression phase . The symbol T* 
denotes the energy associated with the axial motion, whereas T xy is an effec- 
tive transverse temperature. Note that the electrons are quickly heated to 
about 8 keV. 

Similar results are obtained during the current-reduction phase . As shown 
in figure 7, the axial field is now in the opposite direction and accelerates 
electrons toward the center electrode. Figure 11 shows the evolution of the 
electron distribution. The evolution of the electron velocities is shown in 
figure 12. 

The large central region of near-zero magnetic field for the filament model 
allows electron acceleration away from the anode during the compression phase 
and toward the anode during the current-reduction phase. Electron acceleration 
away from the anode as obtained with the filament model is not supported by 
experimental data. The model does explain the experimentally observed electron 
heating. The large number of electrons accelerated to relativistic velocities 
in the region interior to the ring of filaments would give rise to strong mag- 
netic fields. These fields would impede the electrons and lead to a reduction 
in the current. Therefore, a self-consistent computation seems necessary for 
accurate modeling. 


Uniform-Current Model 

The uniform-current model (ref. 14) is probably the simplest model one can 
assume for the plasma focus. Since, in this case, the magnetic field increases 
linearly with r, the region of B = 0 available to produce runaway electrons 
is much smaller than that in the filament model. Thus, the results are much 
less likely to be in conflict with a self-consistent computation. Figure 13 dem- 
onstrates the evolution of the x-y positions of 3000 electrons with an initially 
uniform distribution during the compression phase . The direction of the E x B 
drift is away from the axis, so that the electrons leave a "hole" on the axis. 
Particles are most likely to be accelerated near the axis where B is small; 
but since the induced electric field during the collapse phase drives electrons 
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away from this region, very little acceleration takes place. There is essen- 
tially no heating of the electrons because of trapping by the strong magnetic 
field. The inclusion of collisions would not affect these results appreciably 
since there are at least 10 4 cyclotron periods per collision time. Figure 14 
shows the time evolution of the axial velocity as a function of radial position 
during the compression phase. As can be seen, only a few electrons are heated 
to moderate temperatures near the central region. 

The evolution of the x-y distribution of electrons for the^current- 
reduction phase is shown in figure 15. During this phase the E * B drift 
forces electrons toward the axis. The corresponding evolution of the electron 
velocities as a function of radius is shown in figure 16. Away from the central 
region electrons are effectively tied to the magnetic field and little heating 
or acceleration takes place. However, as electrons forced by the E * B drift 
approach within a Larmor radius of the axis they are free to accelerate under the 
strong eE z force. These electrons quickly gain kinetic energy in the direc- 
tion opposite to J z . Figures 17 and 18 show two examples of such electron tra- 
jectories as they approach the axis and are then accelerated in the z-direction 
toward the center electrode. At a large radius r, these trajectories are well 
described by the drift approximation 

^ . E(r) x g(r) 

B 2 (r) 

However, as the electron approaches within a Larmor radius of the axis, it "runs 
away," although its angular momentum requires it to keep circling the axis as it 
accelerates. The time-averaged energy distribution of the runaway electrons is 
found to be inversely proportional to the electron energy. 

Figure 19 shows a comparison of the X-ray emission experimentally deter- 
mined from the 25 kJ 20 kV plasma focus at the Langley Research Center with the 
emission obtained from the "runaway" electrons impinging on the center elec- 
trode. The X-ray pattern is determined from the electron velocity vectors by 
using the thick target approximation; that is, the emitted energy is uniformly 
distributed in frequency up to the energy of the incident electron. For the low- 
energy emission there is good agreement; whereas, because of insufficient experi- 
mental data, the agreement is questionable for the high-energy emission. 


Runaway-Current Limitation 

In this section the uniform-current model is used to predict some proper- 
ties of the runaway electrons that may be subject to experimental observation. 
For a uniform-current distribution, 

B e (r) = Br 
V n 

ft = — £ en V 
p 2 e v av 

where V OTr is the average drift of the species within radius r and the plasma 
av 

density n g is assumed constant. Particles with Larmor radius greater than 
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their radial position are free to run away. Since the Larmor radius decreases 
with radial position, all particles within a radius r r> determined by setting 
the radial position equal to the Larmor radius, are free to run away. That is, 

r . 1th _ Xth^ 

r ’ eBr r 

where V th = (T/m)^ 2 , the average particle thermal velocity, and is the 

cyclotron frequency. Simplified, the runaway radius becomes 


r 


r “ 



(15) 


If T e = Tj_, the number of runaway electrons will be (m e /m^)^ 2 times the num- 
ber of runaway ions (ref. 13). However, since the ion response time is much 
larger than that of the electrons, ion runaway current would be important only 
on a much larger time scale than that investigated here. For typical focus 
parameters (T = 1 keV; 3=2* 10^ T/m) , the electron r p is approximately 
equal to 10”^ m (see figs. 17 and 18) which is quite small compared to the radius 


of the focus plasma. However, after runaway of electrons occurs within r 


e 


could possibly reach values as high as 3 x 10® T/m (e.g., for n = 10 2 ® m”3 and 
V av = c) resulting in a minimum value for r r of approximately 10 _t> m. 


Further, equation (15) can be used to estimate the total current carried by 
runaway particles: 


I r = neV irr 2 = \/Tm (16) 

av r n 0 e 

This is a remarkable result because it depends upon only a simple experimental 
parameter, the species temperature. For a typical T g = i keV, I r = 400 A 
which is, of course, very small compared to the total current of approximately 
1 MA discharged through a focus. 

The derivation of equation (16) is similar, but not equivalent, to a simple 
derivation of the well-known Alfven-Lawson criterion for the maximum current 
which can be carried by a charge-neutralized, relativistic electron beam 
(ref. 22). In both cases the argument hinges on equating the electron Larmor 
radius to the volume containing the high-speed electrons. The important differ- 
ence is that, in the relativistic-beam case, this velocity is taken as the drift 
speed of the beam V ay ; whereas in the plasma focus, the electron thermal veloc- 
ity V th is the appropriate quantity and V av drops out of the final expres- 
sion for I r . 


DISCUSSION 

Although the computations presented for the filament and uniform-current 
models are not self-consistent, equation (16) is a self-consistent result within 
the assumption of uniform current. That is, if the electron current density 
(for some reason) increases, 6 also increases and r r decreases. Thus, fewer 
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electrons are free to run away and the effects of increasing J and decreasing 
runaway radius compensate each other so that I r remains constant. 

It is important to note that, because the present model assumes that the 
electric and magnetic fields are given, it cannot predict the energy distribu- 
tion of the runaway electrons nor the duration of the associated X-ray pulse . 
Therefore, it is difficult at present to make quantitiative comparisons with 
experimental results. These results are usually given as total runaway electron 
energy inferred from measurements of the hard X-ray energy emitted by the focus. 

Harries et al. (ref. 19) have observed X-rays above 50 keV emitted for 
100 ns and have determined an electron-beam energy of about 1 J. By assuming an 
average electron energy of about 50 keV, they have estimated an average current 
of 200 A, in approximate agreement with the present prediction. However, Gul- 
lickson and Barlett (ref. 18) infer a total energy of 1 kJ for electrons with 
energies greater than about 100 keV in 10 ns, and Maisonnier et al. (ref. 23) 
estimate 600 J for 100 keV electrons. Since these latter results imply runaway- 
current orders of magnitude larger than the present prediction, further work, 
both experimental and theoretical, would be desirable to resolve this 
discrepancy. 

In spite of these limitations, the present model does make an unambiguous 
prediction which is subject to direct experimental verification: that is, run- 

away electrons must be accelerated toward the anode. The reason for this, as 
first discussed in reference 13, and as indicated here, is that axial electric 
fields pointing against the current flow cause charged particles to^drift toward 
larger radii, away from the runaway region near the axes; whereas E parallel 
to 3 "pinches” electrons and ions toward the B = 0 region where runaway can 
occur . 

This prediction implies that ions are accelerated away from the anode, an 
often-demonstrated experimental result, and that electrons run away toward the 
anode. Recent X-ray data (refs. 18 and 19; also, fig. 2) have given strong sup- 
port to this prediction for the electrons. The present prediction is directly 
opposed to that of Newman and Petrosian (ref. 20) who claim that electrons are 
accelerated away from the anode. 

The present work suggests that the X-ray emitting, axially oriented, 
filament-like structures observed by Bostick et al . (ref. 16) are narrow channels 
of high-speed electrons accelerated by electric fields induced during the 
current-reduction phase of the focus. The smallness of r r does not preclude 
the simultaneous existence of more than one of these channels, as suggested by 
Bostick, but the present theory is not yet adequate to predict the conditions 
necessary for the development of multiple filamentation . 

Equations (15) and (16) may have relevance to ions as well, but these equa- 
tions should be applied with care since r r for ions can be of the same order 
of magnitude as the plasma radius, and the uni form- current approximation is prob- 
ably not valid over the whole cross section of the focus. Thus, equation (16) 
suggests, rather than predicts, that runaway ions can carry substantially more 
current and energy than runaway electrons and that the focus may be a good 
device for accelerating very heavy ions. 
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CONCLUDING REMARKS 


Results are presented of a numerical integration of the three-dimensional 
relativistic equations of motion of electrons subject to given electric and mag- 
netic fields deduced from experiments. Fields due to two different models are 
investigated. For the first model, the fields are those due to a circular dis- 
tribution of axial current filaments. As the current filaments collapse toward 
the axis, large azimuthal magnetic and axial electric fields are induced. These 
fields effectively heat the electrons to a temperature of approximately 8 keV 
and accelerate electrons within the radius of the filaments to high axial veloc- 
ities. Similar results are obtained for the current-reduction phase of focus 
formation. For the second model, the fields are those due to a uniform current 
distribution. Both the current-reduction and the compression phases were stud- 
ied . There is little heating or acceleration of electrons during the compres- 
sion phase because the electrons are tied to the magnetic field. However, dur- 
ing the current-reduction phase, electrons near the axis are accelerated toward 
the center electrode and reach energies of 100 keV. The computations are in gen- 
eral agreement with experimental results. A criterion is obtained which limits 
the runaway electron current to about 400 A. The criterion depends only on the 
electron temperature and thus is subject to experimental observation. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
December 15, 1976 
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Figure 1.- Cross-seetional view through coaxial plasma- focus apparatus. 
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(a) Experimental setup. (b) X-ray photograph. 


Figure 2.- Experimental setup to determine direction of electron acceleration in the plasma focus. 
No X-rays are observed from the aluminum target , indicating that electrons are accelerated 
downward toward the center electrode from which strong X-ray emission is observed . 
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Figure 8.- Evolution of the x-y distribution of 3000 electrons during compression phase for the 
10- filament model. Electrons are trapped and heated near and outside the filament radius, 
whereas electrons in the central region are accelerated upward away from the center electrode. 
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Figure 9.- Evolution of the distribution of 3000 electrons in V-r space during the compression pha 
As electrons reach a distance of 1.5 cm from their initial z r 0 axial position, they are 
removed from the calculation. 





Figure 10.- Time variation of transverse temperature T and z-directed 
beam energy T* of electrons during the compression phase in the cur- 
rent filament model . 
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Figure 11.- Evolution of the x-y distribution of 3000 electrons for the current reduction phase 

of the filament model. 
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Figure 12.- Evolution of the distribution of 3000 electrons in V-r space during the current reduction 
phase in the current filament model. Note that electrons with r < R undergo essentially free- 
streaming acceleration with considerable heating for electrons with r near R, The induced E 
is approximately equal to 6 MV/m in the central region. Z 





t = .20 ns 


t = .25 ns 


x, mm 


-4 


-2 


_L 


0 

x, mm 



Figure 13 - - Evolution of the spatial distribution of 3000 electrons during the compression phase for 
a uniform current distribution. The E * 3 drift expels electrons from the central region. 
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Figure 14.- Time evolution of 3000 electrons in V-r space for the uniform current model. Electrons 
are expelled from the central region and are trapped by the strong magnetic field. 
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Figure 15.- Evolution of the spatial distribution of electrons in the current reduction phase of the 
uniform current model. Electrons are forced toward the axis by the E * § drift. 
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Figure 16.- Time evolution of the velocity of 3000 electrons as a function of radius. The 
E x B drift now forces electrons toward the axis. Electrons away from the axis are tied 
to the magnetic field , whereas electrons within r r of the axis are accelerated in the 
negative z-direction. 



(a) Trajectory in x-y space. 


Figure 17.- Typical electron trajectory near the central region during the current reduction phase 
in the uniform current model. Under the influence of the E * B drift, the electron moves 
toward the axis. As the electron approaches within a Larmor radius of the axis, it is acceler- 
ated in the negative z-direction and "runs away." 
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(b) Trajectory in r-z space. 


Figure 18.- Concluded. 
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(a) Low energy (<8 keV). (b) High energy (>15 keV). 


Figure 19.- Comparison of X-ray emission obtained from runaway electrons (solid curve) with 
experimentally determined X-ray emissions (dots) from the 25 kJ 20 kV plasma focus. 
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